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Abstract 

A prototype multi-channel laser Doppler anemometry (LDA) 
processor was assembled using a widc-band transient recorder and a 
microcomputer with an array processor for fast Fourier transform 
(Fni computations. The prototype instrument was used to acquire, 
process, and record signals from a three-component wind tunnel 
LDA system subject to various conditions of noise and flow turbu- 
lence. The recorded data was used to evaluate the effectiveness of 
burst acceptance criteria, processing algorithms, and selection of 
processing parameters such as record length. The recorded signals 
were also used to obtain comparative estimates of signal-to-noise 
ratio between time-domain and frequency-domain signal detection 
schemes. These comparisons indicate that the FFT processing 
scheme allows accurate processing of signals for which the signal- 
to-noise ratio is 10 to 15 dB less than practical with counter 
processors. 

Nomenclature 

D diameter of axisymmetric jet 

f L filter passband low-frequency cutoff 

fn filter passband high-frequency cutoff 

f 0 center frequency of doppler burst 

K amplitude of signal 

L power spectral density of LDA signal 

n noise signal 

N noise apiitude factor 

integrated noise power 
Ps integrated signal power 

(S/N)p frequency-domain signal-to-noise ratio 

(S/Nyp time-domain signal-to-noise ratio 


* Aerospace Engineer. 
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T sampling window time interval 

p LDA burst time constant 


L Introduction 

Laser Doppler anemometry (LDA) offers an indispensable 
method of obtaining nonintrusive, quantitative measurements of 
mean and turbulent velocity characteristics of a broad range of 
Hows, An important measure of the performance of the optical and 
detection systems of the LDA is the signal-to-noise ratio (S/N) of the 
photodetcctor output. This signal must be analyzed to obtain an 
accurate estimate of the burst frequency. Time-domain LDA signal 
counter- processors provide accurate Doppler burst frequency mea- 
surements at high data rates when the S/N exceeds a threshold of 
approximately 10-15 dB. This S/N is often difficult to achieve with 
long focus, backscatter optical configurations, and other factors 
such as surface glare and seeding density variations (due to 
centrifugal flows, for example) compound the difficulty. Since the 
level associated with the white-noise component of the photo-detec- 
tor increases with processor bandwidth, the S/N can be maximized 
by narrow-band filtering, but this may limit the velocity range and 
introduce biasing in flows with severe mean velocity gradients and 
turbulence intensities. Another problem associated with counter- 
processors is discretization or quantization error which masks tur- 
bulence intensities below about 0.5%, depending on the optical and 
electronic set-up. 

Proposals for new processor designs which overcome the 
limitations of countertype processors have been recently made. 
Meyers 1 has proposed a frequency-domain processor which deter- 
mines the burst frequency by interpolating the output of a bank of 
five digital filters. Meyers reports that the processor, which will 
sample at rates of up to 1 GHz, can process bursts with as few as 
150 photons/burst and has a "residual" turbulence level of 0.2% for 
signals for which the burst frequency deviates from the mean by 
20% or less. 

It is also possible to digitize the signal at a high fixed rate 
and execute a discrete Fourier transform to determine the burst fre- 
quency. This approach has been applied by several groups, 2 * 3 - 4 * 5 
and offers several benefits relative to counter-type processors which 
include insensitivity of S/N to processor bandwidth and greater 
noise immunity. Modarress and Tan 4 found that data acceptance 
rates using DFT processing improved by two orders of magnitude 
over counter processors. The method takes advantage of commer- 
cially available transient recorders with sampling rates of up to 
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within the upper 40 dB of the 76 dB dynamic range. A block dia- 
gram of the instrumentation is shown in Fig. 2. The front-end 
detectors of the counter-processors were used to provide the logic 
signals required to initiate digitization of the signal. 

IV. Processing Software 

Test software included modules for acquisition, FFT card 
control, frequency data acceptance/rejection, and velocity component 
computation. A summary of module functions is included in the 
Appendix. In addition to programs for evaluating ensemble aver- 
ages of the flow characteristics, test programs were written to 
evaluate burst S/N ratios for both time domain (counter-processor) 
and frequency-domain (FFT) processing, and to determine the effect 
of sample length on the accuracy of burst frequency measurement. 
This program was run with experimental data to evaluate the effec- 
tiveness of polynomial interpolation of the burst spectra to improve 
estimates of burst frequency. 

A processing rate goal of 25 samples/sec was determined to 
be consistent with minimum wind tunnel sampling periods of 1 min, 
allowing 1500 samples/location. Initial versions of the program 
achieved five accepted sample pairs/sec. Current efforts to optimize 
the code are expected to achieve the sample rate goal without incor- 
porating available hardware accelerators. 

Y,J*ssum 

Fast Fourier transform records for each burst, similar to 
Fig. lb, were used to compare time-domain and frequency-domain 
detection schemes. A value for the time-domain S/N (as seen by 
the counter- processors) was estimated by calculating 10 times the 
common logarithm of the ratio of the sum of burst spectral compo- 
nents to the sum of the nonburst spectral components. The corre- 
sponding frequency-domain S/N value is 10 times the logarithm of 
the ratio of the peak burst spectral component to the peak nonburst 
component. This parameter characterizes the requirement of a soft- 
ware detector to extract the peak spectral component correctly from 
the array of spectral component values. 

Figure 3 shows plots of frequency domain S/N v.s. time- 
domain S/N for two channels with differing noise levels. Typi- 
cally, time-domain S/N values of 0 dB and 10 dB correspond to 
frequency-domain S/N values of 10 dB and 20 dB respectively. 
These values arc somewhat lower than predicted in the previous 
analysis since the noise was integrated over the entire sample length, 
rather than the burst duration. For successful time-domain detec- 
tion, a minimum detectable S/N value of 10-15 dB is required. For 
frequency-domain detection, a comparable difference between the 
peak burst level and the highest noise component is needed. This 
suggests that bursts can be correctly processed in the frequency 
domain which arc 10-15 dB lower than accurately measurable with 
counter processors. Temporary storage of the digitized burst wave- 
forms also permits more rigorous acceptance criteria than those 
employed with the counter processor, allowing greater confidence in 
the accepted data. 

Fast Fourier transform computation time is proportional to N 
log(N), where N is the number of transform points. Efficient sort- 
ing and validation algorithms are comparably scaled. Thus higher 
throughput may be obtained at the expense of frequency resolution. 
The burst data was processed with FFT sizes of 1024, 512, 256, 
128, 64, 32, and 16 points to determine the effect of record length 
on processor accuracy. In addition, the data was also processed 
after setting the first and last 25% of the samples of each burst 
record to zero, simulating the process of padding a shorter time 
series with zeros prior to executing the FFT. 

The peak component of the FFT is an estimate of the actual 
frequency, which may be assumed to be uniformly distributed 
within an interval centered on the measured peak and equal to the 
frequency interval separating the components of the FFT. For this 
assumption, the standard deviation of the estimate is 0.289 times the 
frequency interval, a value which was confirmed by experiment. It 
is possible to improve this estimate by fitting a parabola to the three 
highest spectral levels, and then calculating the peak of the fitted 
parabola. The process is illustrated in Fig. 4, which shows the 
results of a 1024-point FFT and a 16-point FFT with a fitted 
parabola. In each case, the burst envelope was computed and the 
data was extracted from the burst center. The frequency error of 
reduced length estimates was taken as the difference between the 
reduced length estimate and the peak of the 1024-point FFT. The 
effect of sample size and S/N on this error is discussed below and in 
Figs. 5 and 6. 


Figure 5 shows the locus of 50 sample pairs for the 32- 
sample FFT in coordinates of frequency-domain S/N and standard 
deviation of frequency estimate, in which the symbols 1 and 2 rep- 
resent samples from the first and second velocity channels, respec- 
tively. The scatter of the samples in Fig. 5 suggests that frequency 
resolution is more dependent on record length than S/N for the 
selected sample populati on. A further examination of the results for 
the various record length FFTs showed that neither the simple peak- 
component estimate or the parabolic curve-fit estimate approached 
the high-resolution frequency estimate monotonically with record 
length, ruling out any successive approximation scheme. Other 
frequency estimation schemes, such as the maximum entropy 
method- 9 may offer better results but were not considered in this 
study. 

Figure 6 depicts the dependence of the standard deviation of 
reduced length estimates on record length. The three curves corre- 
spond to unpadded peak component estimates, and for parabolic 
interpolations of padded and unpadded records. As can be seen 
from Fig. 6, parabolic interpolation improves the resolution by 
50%.* Replacing the first and last 25% of the data samples with 
zeros had no significant effect on the frequency estimate. This sug- 
gests that burst-length velocity biasing, as described in the back- 
ground section, may be alleviated by zeroing a fixed number of 
samples outside the burst center without adversely affecting fre- 
quency resolution. 

Velocity profiles of a 1.5-in.-diam axisymmetric jet with an 
nozzle speed of 32 m/sec were obtained to test the measurement and 
survey programs, and to detect serious biasing or defects in the 
hardware and software. A copy of a typical data sheet, with his- 
tograms for one location is presented in Fig. 7. Measurements of 
the mean velocity at a station four jet diameters downstream of the 
nozzle exit are presented in Fig. 8, with measurements from a hot- 
wire anemometer and LDA counter-processor for comparison. 
Similar comparisons of the normal and shear stresses for the axial 
and radial components are shown in Fig. 9. These results show 
generally good agreement between the two methods. Improved 
sampling rates will permit larger sample populations and greater 
confidence in the estimation and acceptance algorithms. 

.YL Confiding Remarks 

The present study demonstrated the feasibility of a micro- 
computer-based multichannel LDA signal processor using transient 
recorders and FFT. The throughput of this prototype processor is 
sufficient for most low-speed wind tunnel applications and the per- 
formance is superior to counter-processors for signals with S/N less 
than 10 dB. Preliminary evaluations of several processing algo- 
rithms were based on Doppler burst waveforms obtained from a 
multicomponent LDA. Results of these evaluations suggest that 
zeroing a fixed number waveform samples away from the burst 
center may minimize burst duration velocity biasing without 
adversely affecting frequency resolution, and that parabolic interpo- 
lation of the spectral components may improve the frequency reso- 
lution by about 50% relative to simply choosing the peak compo- 
nent. Continued trends of increasing computational and cost per- 
formance of microcomputer-based instrumentation should improve 
the availability of this approach as a cost-effective alternative to 
dedicated, nonprogrammable LDA signal processors. 
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APP ENP1X ; SQETWARE..QRGANIZA TION 

(software modules in parenthesis) 

A. Initialization 

L Initialize transient recorder and array processor 
(Set_up3, initJTtsys) 

2. Input test parameters: digitizing rate, signal bandwidth 

mix frequency, number of samples 
acceptance criteria 

3. Initialize graphics display (graph Jnit, text Jnit) 

B. Acquisition Loop 

1. Transfer data from transient recorder to computer 
upon receiving event pulse. (Trans2k) 

2. Convert and scale data to 16-bit integer format (Trans2k) 

3. Execute FFT and PSD (rcaLfft_psd) 

4. Son frequency data for largest element (sortl,sort2) 

5. Acccpt/rcjcct on basis of relative magnitude and location 

of next highest elements (reject if non-adjacent) (sortl,sort2) 

6. Update statistical quantities and graphics display (update) 

7. Repeat for each channel until sample number requirement 
is fulfilled 

C. Data processing/display 

1 . Compute turbulence quantities (avg) 

2. Update graphical displays (print_avg) 

3. Store data and move to next traverse location (copy) 
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Figure L (la) LDA single burst waveform, (lb) Single burst PSD. 



Figure 2. Instrumentation diagram (1 of 2 channels). 



Figure 3. Typical variation of (S/N)p vs (S/N>f. 
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Figure 4. Parab olic in terpolation of 16-point FFT and comparison 
with 1024- point FFT. 
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Figure 5. Typical distribution of interpolated frequency error. 
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Figure 6. Effect of sample size on interpolated frequency error. 


Figure 8. Mean velocity measurements of axisymmetric jet 
Uc = 32 m/sec, D = 3.81 cm, x/D = 4.0. 
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Figure 7. Typical output page for two-channel set-up. 



Figure 9. Turbulence measurements of axisymmetric jet. 
Ue = 32 m/sec, D = 3.81 cm, x/D = 4.0. 
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